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ABSTRACT 

The flux variability of blazars at very high energies does not have a clear origin. Flux variations on time scales down to the minute 
suggest that variability originates in the jet, where a relativistic boost can shorten the observed time scale, while the linear relation 
between the flux and its RMS or the skewness of the flux distribution suggests that the variability stems from multiplicative processes, 
which are associated in some models with the accretion disk. We study the RMS-flux relation and emphasize its link to Pareto 
distributions, characterized by a power-law probability density function. Such distributions are naturally generated within a minijets- 
in-a-jet statistical model, in which boosted emitting regions are isotropically oriented within the bulk relativistic flow of a jet. We 
prove that, within this model, the flux of a single minijet is proportional to its RMS. This relation still holds when considering a large 
number of emitting regions, for which the distribution of the total flux is skewed and could be interpreted as being log-normal. The 
minijets-in-a-jet statistical model reconciles the fast variations and the statistical properties of the flux of blazars at very high energies. 

Key words. Relativistic processes - Galaxies: jets - Galaxies: active - Gamma rays: galaxies - X-rays: bursts 
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1. Introduction 

Active galactic nuclei (AGN) are astrophysical sources emitting 
broad-band electromagnetic radiation that can be observed from 
radio wavelengths to very high energies (VHE, E > 100 GeV). 
The current framework for understanding the emission from 
these objects depicts them as composed of a disk, feeding a super 
massive black hole (SMBH), possibly with a jet on each side of 
the accretion plane. The orientation of the system with regard to 
the observer's line of sight is one of th e key parameters for uni - 
fying the various subclasses of AGN jUrrv & Padovanilll995h . 
The class of blazars for which the jet is probably directed within 
a few degrees toward the observer is of particular interest for the 
study of non-thermal processes occurring in astrophysical rela- 
tivistic flow. 

The VHE emission of blazars, such as the BL Lac ob- 
ject Mrk 421 or the flat spectrum radio quasar PKS 1222-1-21, 
can s ometimes exhibit doubling times of the order of ten min- 
utes jGaidos et al.lll996l: lAleksic et ai]|2011h while variability 
on the minute scale is observed in light curves of BL Lac ob- 
jects, such as Mrk 501 d^bertetal. 2007) or PKS 2155-304 
dAharonian et alj|2007h . The high quality of the VHE data on 
the latter object has enabled the characterization of various sta- 
tistical properties of the emission, such as a skewed flux distri- 
bution, interpreted as a log-normal distribution, and a linear re- 
lation between the sample RMS of the flux and th e sample flux 
dH.E.S.S. CoUaboration. Abramowski et alj|2010l) . 

These properties were initially studied in the X-ray 
emission of the galacti c black hole binary Cygnus X-1 
dUtdev & McHardvlboOll: lUttley et al.ll2005h and are now ob- 
served in other accreting o bjects, such as the blazar BL Lac 
dGiebels & Degrangd l2009l) . or n on-aUgned AGNs, suc h as 
the Seyfert galaxi es NGC 405 1 dMcHardv et alj |2004|) and 
IRAS 13224-3809 dGaskellll2004 . 

They are sometimes interpreted as arising from multiplica- 
tive processes, originating e.g. in the accretion disk. The prod- 



uct of random variables whose logarithms have finite moments 
is indeed asymptotically log-normal. Furthermore, flux log- 
normality implies a linear relation betw een the flux and its RMS 
(see Appendix D of lUtdev et ani2005h . 

The statistical properties of non-thermal light curves are key 
ingredients for answering the crucial question of the origin of 
variability in AGNs: do the variations of the flux come from the 
jet itself, i.e. the medium that conveys the accelerated particles, 
or from the disk that could modula te the jet emission . In ac- 
cretion models such as developed by iLvubarskiil dl997h . inward 
propagating fluctuations in the d isk from diff'erent radii c umulate 
in a multiplicative way (see e.g. lArevalo & Uttley||2006l ). which 
could explain the log-normality and the linear RMS-flux rela- 
tion. 

However, variability on time scales shorter than the light- 
crossing time Tr, ^ IGMjc^ are difficult to achieve 
dNaravan & Piranl l2012h in disks. Transient emission on time 
scales fvar ~ TqI50 are easier to locate in a jet with a bulk 
Lorent z factor F (unless gene rated by a companion system, as 
e.g. in iRieger & Volpd l2010h . Phenomena occurring on time 
scales Tg then appear at fvar ~ Tc/ 6 to the observer, where 
6 is the Doppler factor of the jefl. 

Values of the Doppler factor 6 > 50 are not uncommon 
for VHE blazars and p ermit both minute-time-scale variations 



dAharonian et alJ 2007[) and an optical thinness to the emit- 



ted VH E Y-rays dBegelman et alJl2008l) . lGhiseUini & Tavecchiol 
d2008h propose a needle-in-a-jet model, where fast minijets 
close ly aligned with the li ne of sight could explain the varia- 
tions. O^nniosetaL l (l2009h have investigated a jets-in-a-jet sce- 
nario, involving r econn ection and a Poynting-dominated flow. 
iNarayan & PiranI d2012h considered a turbulent process, mod- 
elled by a random motion of minijets inside the jet. Models 
involving minijets have also been invoked to e xplain the fast 
variations in the emission of radio galaxies dGiannios et al.l 
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' The link between the Doppler factor and the bulk Lorentz factor is 
discussed extensively in the following. 
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120 lOh . though the time scales of variations observed at VHE 
are m uch longer than for blazars, approximately a day for M 87 
jAbramowski et al. 201^. 

These models are additive (the flux is given by the sum of 
the contributions of several regions) and, provided there are 
finite moments for the individual components, the central limit 
theorem (CLT) can be applied. But the Gaussian flux distribu- 
tion expected from the CLT is not skewed, as is the observed one. 

We seek a way out the jet-or-disk origin dilemma by study- 
ing the implications of the observed linearity between the sam- 
ple RMS and the flux (Sect. 12.11 ). We discuss the properties of 
Pareto-distributed fluxes in Sect. 12.21 and discuss the relevance 
of the CLT for these random variables. We show in Sect. [3] that 
Pareto distributions are naturally generated by the minijets-in-a- 
jet statistical model and discuss the consequences of this model 
in Sect.H 



2. Additive processes and Pareto distributions 

Modelling of VHE blazar variability should reproduce the short 
time scale of the variations and the statistical properties of the 
flux. While the stochastic nature of the flux is understood if the 
emission is modulated by the disk of the SMBH, the observed 
time scales are difficult to explain. Inversely, a jet origin of the 
variability is often rejected because the models are additive and 
should not be able to reproduce the skewness of the flux nor the 
linear RMS-flux relation. This argument is studied in the next 
section. 



2.1. Additive and multiplicative processes 

Within additive models, the observed flux is the sum of the 
contributions of several (and potentially many) regions. If the 
components are modelled as independent, identically distributed 
random variables (hereafter iid) with finite moments, their sum 
should follow a normal distribution according to the CLT, as- 
suming that the number of regions is large enough - typically 
more than a few tens. A natural outcome of such a process would 
then be symmetric flux distributions, in disagreement with the 
high-flux tails frequently observed in light curves. Moreover, the 
addition of incoherent variations in various components should 
result in a flux being independent of its variance. A translation 
into statistical terms would be that the sample mear£| and the 
sample RMS of a distribution are independent if and only if the 
underlying dist ribution is a Gaus sian (see the discussion in the 
Appendix D of lUttlev et al1l2005h . 

A contmrio, multiplicative processes naturally generate log- 
normality, with a characteristic tailed and skewed distribution. 
If we let the quantity ^ be the product of a large number of 
iid variables, then the logarithm of (/> is a sum of iid quantities 
and, according to the CLT, log (p follows a normal distribution, 
i.e . follows a log-normal distribution (see e.g. the simulations 
of lloka & Nakamtirall2002l with only three multiplicative vari- 
ables). 

The RMS -flux relation can be expl ained as a consequ ence 
of a log-normal distribution of the flux (lUttlev et al.ll2005h . but 
the reciprocal is not true. Indeed, let the observed flux be a 
function of an underlying random variable x, so that (p - fix)', 
for a log-normal distribution, / is the exponential function and 
X is normally distributed. A small fluctuation of x around xq, 6x, 



results in a small fluctuation of the flux (p around /(xq), 
the variance of the flux is 



3/ 
dx 
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With a sample flux proportional to /(xq), a linear flux-RMS 
relation is equivalent to 
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Equation (|2]i is one of the definitions of the exponential func- 
tion, so the sample RMS is proportional to the sample flux if 
and only if the flux is the exponential of an underlying variable, 
as is the case for a log-normally distributed flux. Log-normality 
thus implies linearity between the RMS and the flu5(0, as does 
any distribution that arises from the exponential of an underly- 
ing variable. 

2.2. Pareto distributions and the central limit theorem 

We now consider the class of distributions that, instead of aris- 
ing from the exponential of normal distributions, arise from the 
exponential of exponential distributions exp(-Q'x). These are 
known as Pareto distributions and are characterized by a proba- 
bility density function (PDF) that follows a power-law function 
of index I + a: 



friy) 



r 



for y > 1 . 



(3) 



Pareto distributions are used in seismology, in finance, and in 
biology and various examples of applications such as "volcanic 
eruptions, solar-flares, lightning strikes, river networks, forest 
fires, extinctions of biological species, war casualties, Internet 



traffic, stock retur ns, insurance pay-off's", can be foun d in the 
literature (see, e.g.. lNewrna nll2005HZaliapin et aUlIOOl and ref- 
erence therein). 

The CLT cannot be applied to Pareto distributed iid, mostly 
because, for a < 2, the variance of each random variable cannot 
be defined. A generalized central limit theorem has been estab- 
lished for such a high-tailed distribution (see, e.g., IVoilll2005l) . 
and it states that the sum of Pareto-distributed iid converge to a- 
stable distributions, a class of probabiUty distribution with fasci- 
nating properties. 

Only a handful of c-stable distributions can be expressed 
in terms of elementary functions. Some special cases are sym- 
metric, such as the Gaussian distribution, which corresponds to 
the limit a - 2, but a large skewness is generally seen, e.g., 
with the L andau distribution (a = 1), used by particle physicists 
dGroom & Klein 2000) to describe the distribution of the energy 
loss of a charged particle passing through matter. In particular, 
for < a < 2 the PDF of an g-stable distribution asymptotically 
follow s a power law of index 1 + a (IZaliapin et al.ll2005t iNolanI 
I2012h . 

The convergence of the CLT is achieved with a few tens of 
"regular'Q components, but sums of Pareto iid variables tend to 
the asymptotic or-stable di stributions only for l arge numbers of 
components, typically 10"^ (IZaliapin et al.ll20()5l) . which prevents 
the PDF of the sum from being analytically derived for an in- 
termediate number of Paretian iid quantities. Even though the 
log-normal distribution is not an attractor for the sum of such 



i.e. averaged on short time periods in the light curve. 



And the reciprocal is not true! 
i.e. with finite moments. 
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iid variables, skewed distributions that closely resemble the log- 
normal can be expected, as is shown in Sect. 13.31 

Equations ([TJ and ^ prove that the sample RMS is propor- 
tional to the sample flux for a single random variable as long as 
it is the exponential of an underlying variable. This is the case 
for Pareto distributions, and it will be shown in Sect. l3.3l that this 
relation holds for the sum of a large number of Paretian iid, cer- 
tainly because the PDF of the sum follows the same power-law 
behaviour for high fluxes. 

3. A minijets-in-a-jet statistical model 

Power laws are generally seen as the signature of a scale invari- 
ant process and a wide variety of models can certainly produce 
such flux distributions. In the next sections, a simple kinematic 
model that naturally generates Paretian fluxes is developed. 

The relativistic enhancement of the flux, known as Doppler 
boosting, is computed for an emitting region randomly oriented 
in a medium. For the sake of generality, both the emitting re- 
gion and the medium, called minijet and jet respectively, have 
a relativistic motion. Beamed astrophysical objects thus fit in 
this general framework, as well as non-beamed objects, when 
the Lorentz factor of the medium is set to unity. 

3.1. Doppler factor of a minijet randomly oriented in a jet 

We consider a minijet, e.g. a blob of plasma, that moves ran- 
domly in a jet and is characterized by its Lorentz factor y and as- 
sociated velocity /3 - ^[T^-Y^^ within the jet frame. The orien- 
tation of the minijet in the jet frame can be defined with the con- 
ventional angles in spherical coordinates, which will be called ijj 
and If. The jet is defined by a Lorentz factor F, associated veloc- 
ity E, within the observer frame and its axis is misaligned from 
the line of sight by an angle 9. A schematic representation of 
these parameters is shown in Fig.[T] 

To derive the Doppler factor 6 of the emitting region, the 
minijet, one only needs to determine the energy measured in 
the observer frame as a function of the energy in the minijet 
frame, where the emission is assumed to be isotropic. We con- 
sider a photon emitted in the minijet frame with energy iSminijet 

and momentum Pminijet = [P.x minijet, Py minijet, Pz minijet}- A SCricS § ; 

of transformations must be applied to derive the energy in the 
observer frame : the tilt of the jet, a boost, the tilt of the mini- 
jet in the jet, defined by two angles, and finally the boost of 
the minijet. Only photons propagating along the Xobs direction, 
i.e. along the line of sight, are considered because these are the 
only ones observed. We can thus set the transverse momentum 
to zero, Py obs = Pz obs = 0, while the longitudinal momentum is 
Px obs = E, where the speed of light is set to unity and where E 
is the photon energy in the observer frame. 

These equalities and the frame transformations are summa- 
rized in Eq.(|4]i, where c and s denotes the cosine and sine func- 
tions: 
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The inverse of the Doppler factor, i.e. the ratio between 
the emitted and the observed energies, is derived by inverting 
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Fig. 1. Schematic view of the geometry. The left side corre- 
sponds to the observer frame ^^obs, in which the jet is tilted by 
angle from the line of sight (along Xobs) and is boosted by a 
Lorentz factor F (velocity E) along x' . The minijet is defined 

jet 

by its Lorentz factor y (velocity /?), in the jet frame ^jet (top 
right). The orientation of the minijet along xV^ijgj compared to 
the jet axis in its rest frame, xjet, is defined by the spherical an- 
gles tp and (/?. The emission is assumed isotropic in the minijet 
frame 'Rminijet (bottom right). 



Eq.(|4]i. The time-like component of the vectorial equality then 



reads as 6 E 



^minijet' 



1 



with 6 given in Eq.(|5]i: 



rF(l + Y^pc^ - (E +y6c^)Cfl) + ypses^c^ 



(5) 



For a jet closely aligned with the line of sight, the extrema of 
the Doppler factor are 



6< 
6 > 



yr(i-2;)(i-/j) 
1 



4ry 



(6) 



where the approximated expressions correspond to the ultra- 
relativistic limit. 

It is worth noting that for reasonable Lorentz factors such as 
F = 5 and y - 5, the maximum Doppler factor obtained for a 
minijet and a jet aligned with the line of sight is enhanced by a 
factor 2y = 10 compared to isotropic emission in the jet frame 
(y = 1), and Doppler factors as high as 10^ are reached. 

Figure |2] shows the Doppler factor as a function of the mis- 
alignment from the line of sight 9. The Lorentz factors of the jet 
and of the minijet are fixed to a value of five and 9 is normalized 
to 1 /F, the opening angle of the jet. 

The maximum and minimum values of the Doppler factor 
are represented with the long and short dashed lines. The average 
Doppler factor is derived with fi = cos tfr uniformly distributed 
between [-1,1] and ip uniformly distributed between [0,27r]. This 
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lO"' 1 

Misalignement from the line of sight 6 / 

Fig. 2. Doppler factor of a minijet with a Lorentz factor y - 5 
moving in a jet of Lorentz factor F = 5 as a function of the mis- 
alignment between the jet and the line of sight 9, normalized to 
the opening angle of the jet 1 /F. The maximum, average, and 
minimum Doppler factor are compared to the usual Doppler fac- 
tor (solid line), derived for an isotropic emission in the jet frame 

(r= !)■ 

corresponds to an isotropic distribution of the minijet orientation 
in the jet frame. One can compute the average Doppler factor in 
the blazar case, i.e for 6 = 



(5o = 



6o = 



rF(l 



1 



1 



In 



rF(i-s)^ \i 



— In 2y. 

y 



(7) 



(8) 



For y > 4, the reference Doppler factor (y = 1), represented 
with the solid line in Fig. |2] is larger than the average Doppler 
factor, a comparison that can be extended only to quantities that 
are linearly dependent on the Doppler factor. The impact on the 
flux intensity, which is roughly a quartic function of the Doppler 
factor, is discussed in the following. 

3.2. Flux of a single minijet 

Let I{E) be the flux intensi ty at energy E, the quantity 1(E) /E^ 
is a Lorentz invariant (e.g.. lRvbicki & Lightmanl[l979h . and the 
Doppler boost 6 of the energy results in an enhancement 6^ o f 
the intensity or of the flux (see, e.g., lUrrv & Padovanilll995l) . 
One finds the usual enhancement of the bolometric luminosity 
by a factor d'^ after integration over the energy domain. 

The non-thermal spectra of high-energy astrophysical 
sources can generally be approximated, at least locally, by 
power-law functions. We characterize the intensity of the source 
by the index s (photon index s + 1), i.e. /(fiminijet) « E^.^..^^, 
where /iminijet is the photon energy in the source rest frame. The 
intensity measured by the observer is, in the observer's frame. 



1(E) = 5'/(£™nijet) = 6'I(E/5) <x 6'(E/6)-' ^ 6^^'E 



(9) 



The distribution of the intensity can then be derived from 
Eq.(|5]l, imposing a distribution for the underlying random vari- 
ables. We assume in this model that the relevant physical param- 
eter is the orientation of the minijet(s) in the jet frame. Imposing 



isotropy in this frame corresponds to having if/ and ip angles ho- 
mogeneously populate the unit sphere, i.e. yU = cos i/' uniformly 
distributed between [-1,1] and (p uniformly distributed between 
[0,2;r]. For the sake of simplicity, we derive the probability den- 
sity function of the intensity in the blazar case. In this case, com- 
bining the definition of the Doppler factor in Eq.(|5]) with Eq.(|9l) 
yields 



1(E) cx [yF(l - Z)(l -/3i.i)]-'-'E-' = (4Ty)'^'g(M)E-\ (10) 
The factor (4F7')^^* scales the function g so that 

< 1 withyL/e [-1,1]. (II) 



{(l+X)(l+^ 1-/3 
g(M) = I : X 



We call In = g(fj) the intensity normalized to its maximum. 
The PDF of In, fAh), is linked to the PDF of // = cos i/r, fc(n) = 
1/2 with e [-1; 1] via 



Mn) = 



dg-\lN) 



dl 



fc[g'\lN)), 



(12) 



1 / 1+Z 1 __L 

where g \x) = ^ x 1 1 — ^ ^ ^ ''^ 

Then, the probability density function of the normalized in- 
tensity is 



Mh) 



1 +E 
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1 



2/3 Ay\3 + s) ,i 



(13) 



The flux emitted by a minijet randomly oriented in a jet 
(or in a non-relativistic medium) thus follows a Pareto distri- 
bution of index a = 1/(3-1-5), where s is the spectral in- 
dex of the intensity. The Pareto distribution holds for a non- 
beamed obje ct (F = 1 and S - 0), as derived independently 
bv lClausen-Brown & Lvutikovl(l2012h to model the flares of the 
Crab in the high-energy domain. 

The Pareto distribution also holds for beamed objects mis- 
aligned from the line of sight, as shown in Fig.|3] where the dis- 
tribution of the logarithm of the intensity oc 5^^'^ has been gener- 
ated by drawing random minijet orientations. The conservation 
of the slope of the distribution can be understood by neglecting 
the term y/^ses^c^ in Eq.( [Sj, in which case the previous proof 
holds with the inverse of the Doppler factor remaining a linear 
function of cos tfr. The conclusions derived in the following can 
then be applied to both blazars, such as PKS 2155-304, and non- 
aligned objects, such as M 87. 



3.3. Fiux of N minijets 



As discussed in Sect. 12.21 Pareto-distributed random variables 
do not follow the central limit theorem, and their sum does not 
converge to a normal distribution. Assuming that the jet is com- 
posed of several randomly oriented minijets, the total emitted 
flux is proportional to the sum of the independent minijet in- 
tensities, in an optically thin medium. 

To illustrate this point, minijets of Lorentz factor y = 5 are 
simulated for = and a jet boost F = 5. These moderated val- 
ues yield a single component flux spanning over eight decades, 
roughly between (F/y)"*^* and (4Fy)-'+'. Unless strongly sup- 
pressing the beaming (typically for y < 1.5), the dynamic range 
of the flux is large enough that the location of the cut-off' does 
not affect the following results (up to at least 10"^ additive com- 
ponents). 
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Fig. 3. Distribution of the logarithm of the intensity of a mini- 
jet for various angles 9 between the line of sight and the jet 
axis. The intensity is proportional to 6^'^", where 6 is the Doppler 
factor and s the spectral index. These distributions are obtained 
simulating random orientations of a minijet for = 1, for a jet 
boost r = 5, and a minijet boost y - 5, without generality loss. 
While the intensity decreases when the misalignment increases, 
the slope of the distribution remains the same as in the "blazar" 
case, 9 = 0. 

As shown in Sect. 13.21 the results for other orientations or 
other Lorentz factors remain identical within a multiplicative 
factor (i.e. a shift in log representation). The flux distributions 
simulated for several numbers of minijets are shown in Fig.H] 
To generate smooth distributions, 10^ iterations, called here- 
after time steps, are performed with 6 {1,10,30,10^,3 x 
10", 10^3x 10^10"*}. 




Log^^ l(N blobs) 

Fig. 4. Distribution of the logarithm of the intensity of inde- 
pendent and randomly oriented minijets. The number of mini- 
jets increases from left to right with A^ 6 {1, 10, 30, 10", 3 x 
10^, 10^*, 3 X 10^, lO'*}. Even for a large number of regions, asym- 
metrical, tailed distributions are obtained. 

For a large number of minijets, typically A^ > lO"*, the dis- 
tribution of the logarithm of the flux can be described with a 
peak, followed by a power-law tail. The distributions obtained 
for 3 X 10^ and 10^ minijets are shown Fig. |5] with a linear y- 
axis, together with the best-fit functions corresponding to a nor- 
mal and log-normal flux. 

Although the distribution of the flux of A^ minijets is neither 
normal nor log-normal. Fig. |5] illustrates how an experimental 



5f 




Log^^ l(N blobs) 



Fig. 5. Distribution of the logarithm of the flux of A^ minijets 
for A^ = 3 X 10^ (left) and A^ = lO'* (right), as in Fig. H The 
continuous black and grey dashed lines represent the best fit with 
a log-normal and normal flux distributions, respectively. 

distribution, with limited statistics and a poor characterization 
of the high flux tail, could be interpreted as a log-normal distri- 
bution, even if arising from an additive process and not from a 
multiplicative one. 

We have shown in Sect. 12.21 that the RMS is strictly pro- 
portional to the flux if and only if the flux is the exponential 
of an underlying variable, which is the case for Pareto distribu- 
tions. Interestingly, linearity holds when considering the sum of 
a large number of power-law components. To illustrate this state- 
ment, light curves are simulated with lO*" time steps for A^ = 1 
and A^ = 10'^ minijets. The sample mean and the sample RMS 
are then computed in ten points wide windows. For the sake of 
clarity, the flux and the RMS are averaged in 50 bins between 
the maximum flux and the minimum flux. The sample RMS is 
plotted as function of the sample flux in Fig. |6] where the error 
bars correspond to the standard deviation in each bin. 

The positive x-intercept in the RMS-flux relation, as seen in 
the right-hand panel in Fig. |6] corresponds to the peak in the 
flux distribution shown in Fig. |4] While it could be interpreted 
as a constant component, this threshold flux corresponds to the 
value around which the average emissions of the minijets pile up 
within the additive model presented here. 

4. Discussion and conclusion 

The linear relation between the RMS and the flux of individ- 
ual blazars or the skewness of the distribution of the flux are 
interpreted by several authors as arising from multiplicative pro- 
cesses, favouring a variability stemming from the disk. 

We are studying a minijets-in-a-jet statistical model where 
the variability stems from the jet itself. We first considered the 
enhancement of the flux due to the relativistic Doppler effect 
and show that an isotropic orientation of a single region in the 
jet frame result in a Pareto distribution of its contribution to the 
flux, characterized by a power-law PDF. As for log-normal vari- 
ables, Pareto-distributed variables can be seen as the exponential 
of an underlying variable and, as such, satisfy the linear relation 
between the flux and its RMS. These random variables do not 
fulfil the hypotheses under which the central limit theorem can 
be applied and the summation of the contributions of many re- 
gions does not yield a normal distribution of the flux. The total 
flux distribution thus remains highly skewed and tailed when in- 
creasing the number of contributing regions, up to at least 10^ 
regions, and could very well be interpreted as a log-normal. A 
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Fig. 6. Sample RMS as a function of the sample flux of = 1 
minijet (top) and the sum of N - 10'^ minijets (bottom). Linear 
relations are found in both cases, with a zero x-intercept in the 
first case and a positive one in the second. 



noticeable difl^erence between the attractor of the sum of Pareto 
distributed variables, called a-stable distributions, and the log- 
normal distribution is their high flux tails that follow a power law 
and that could be probed with long-term, finely-sampled light 
curves. This power-law tail is, moreover, responsible for the in- 
variance of the linear relation between the RMS and the flux 
when the number of minijets is increased. The study of the link 
between the slope of the RMS-flux relation and the kinematic 
parameters of the mini-jets-in-a-jet statistical model is beyond 
the scope of this paper and will be the object of further studies. 

The minijets-in-a-jet st atistical model is lar gely inspired 
from the jets-in-jet model of iGiannios et al. I (l2009h . where mag- 
netic dissipation triggers the emission of relativistic blobs of 
plasma. The authors reproduce the luminosity and the fast vari- 
ations observed in the light curves of blazars at very high en- 
ergies with a Lorentz factor of the minijet y ~ 10 in the jet 
frame. iNaravan & PiranI (1201 2h refined this model and invoke 
relativistic turbulence to explain both the variability time scales 
and the scarcity of the flaring events. The random orientation of 
the minijet in the jet frame would be directly linked to the direc- 
tion in which the reconnection region is created (magnetic dis- 
sipation) or to the wandering of the velocity vector of the mini- 
jet over the radiating time (relativistic turbulence). The isotropic 
orientation of the minijets in the jet frame is the most natural 
configuration in both scenarios. 



A " minijet statistical model'B has als o recently been de- 
rived bv lClausen-Brown & Lvutikovl(l2012h to explain the high- 
energy flares of Crab Nebula. The authors assume that magnetic 
reconnection triggers the emission of two-sided relativistic flows 
with a random orientation in the nebula frame, a particular ge- 
ometry corresponding to a bulk Lorentz factor of the medium set 
to unity in the framework developed in this paper Assuming that 
the individual components do not overlap, the authors restricted 
the study to a very small number of minijets, a problem that is 
overcome in this paper. 

Unlike Iciausen-Brown & Lvutikovl(l2012h . we did not study 
observables such as the power spectral density or the variation 
time scales observed in light curves. These quantities depend 
on the detailed mechanism responsible for the generation of the 
minijets, and their study is beyond the scope of this paper. Our 
model, based on a Doppler-boosted beamed emission, can cer- 
tainly be generalized to other mechanism s, such as the "kinetic 
beaming" studied by ICerutti et alJ (1201 2h . which will probably 
yield similar conclusions on the statistical properties of the flux, 
while providing refined energy-dependent observables. 

The generalization of such stochastic models, where many 
sub-regions are considered, growingly reduces the dominance of 
shock emission models in high-energy astrophysics. The spring- 
board of such a change of paradigm points out the need for high- 
quality non-thermal light curves and widening studies of flux 
distributions, of the relation between the sample RMS and the 
sample flux, and of the power spectral density. 
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